Stability of periodic traveling waves for the 
quadratic and cubic nonlinear Schrodinger equations 



Sevdzhan Hakkaev^, Iliya D. Iliev^ and Kiril Kirchev^ 
December 20, 2011 

^Faculty of Mathematics and Informatics, 
Shumen University, 9712 Shumen, Bulgaria 

^Institute of Mathematics and Informatics, 
Bulgarian Academy of Sciences, 1113 Sofia, Bulgaria 

Abstract 

We study the existence and stability of periodic traveling-wave solutions for 
the quadratic and cubic nonlinear Schrodinger equations in one space dimen- 
sion. 

1 Introduction 

In this work we consider the nonlinear Schrodinger equation (NLS) 

iut + Uxx + \u\^u = 0. (1.1) 

This equation appears in various problems, modeling many phenomena such as the 
behavior of a non-ideal Bose gas with a weak particle interaction, the spreading of 
the heat impulse in solids, the Langmuir waves in a plasma, etc. [171 ITS] . 

Our principal aim is to study the the orbital stability of the family of periodic 
traveling-wave solutions 



u 



ifi{x, t) = e*(^^+("^-''')*V(x - 2vt). (1.2) 



where r{y) is a real-valued T-periodic function and f,^ G M are parameters, for 
quadratic (p=l) and cubic (p=2) nonlinear Schrodinger equation. The problem of 
the stability of solitary waves for nonlinear dispersive equations goes back to the 
works of Benjamin [4| and Bona [5] (see also [H [151 IIS])- A general approach for 
investigating the stability of solitary waves for nonlinear equations having a group of 
symmetries was proposed in [8] . The existence and stability of solitary wave solutions 
for equation (11.11) has been studied in [12]. Recently in [3], the authors developed a 
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complete theory on the stabihty of cnoidal waves for the KdV equation. Other new 
exphcit formulae for the periodic traveling waves based on the Jacobi elliptic func- 
tions, together with their stability, have been obtained in [21 El [ID] for the nonlinear 
Schrodinger equation, modified KdV equation, complex modified KdV equation, and 
generalized BBM equation. In [11], the stability of periodic traveling wave solutions 
of BBM equation which wave profile stays close to the constant state u = {c— 1)^/^ 
is considered. 

In this paper, we prove stability of the periodic traveling waves (11. 2p not oscillating 
around zero (r 7^ 0) for the quadratic and the cubic NLS. Our main results are 
formulated as Theorem 12.11 and Theorem 13.11 in Sections 2 and 3 below. We base our 
analysis on some appropriate invariant laws. Our approach is to verify that (f is a 
minimizer of a properly chosen functional M which is conservative with respect to 
time over the solutions of fll.ll) . We consider the L^-space of T-periodic functions in 
rc e M, with a norm ||.|| and a scalar product (.,.). To establish that the orbit 

= {eX._^,t): (e,r;)e[0,T]x [0,277]} 

is stable, we take 

t) = e'^'ifix - ^, t) + h{x, t) = e^^ [r{x - ^ - 2vt) + h + i/is] 

and express the leading term of M{u) — M{ip) as {Lihi, hi) + (-^2/12, ^2) where Lj are 
second-order selfadjoint differential operators in L^[0,T] with potentials depending 
on r and satisfying Lir' = = 0. The proof of orbital stability requires that zero 
is the second eigenvalue of Li and the first one of L2. 

Recall that the quadratic and cubic nonlinear Schrodinger equations are globally 
well-posed in i7*(T), for s > [BJ. 

The paper is organized as follows. We consider the quadratic and the cubic cases 
in Sections 2 and 3 respectively. In Appendices 1 and 2, some facts from the theory of 
complete Abelian integrals (e.g. Picard-Fuchs equations, polynomial moduli) are used 
in order to derive several inequalities we needed during the proof our main results. 



2 Existence and stability of periodic traveling waves 
for the quadratic Schrodinger equation 

Consider the equation 

iut + Uxx + \u\u = 0, (2.1) 

where u is a complex-valued function. 

We are looking for a solution of equation (12.11) in the form (ll.2p where r is real- 
valued. For r one obtains the equation 

r" -ur + r\r\ = 0. (2.2) 

Therefore, 

2 

r — oor -\ — r \r\ = c (2.3) 
3 
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and r is periodic provided that the level set H{x, y) = c of the Hamiltonian system 
dH = 0, with 

2 

H(x, y) = — ojx^ + 

contains a periodic trajectory (an oval). The level set H{x,y) = c contains two 
periodic trajectories if a; > 0, c G (— |c(;^,0) and a unique periodic trajectory if 
w G M, c > 0. Under these conditions, equation (12.31) becomes H{r,r') = c and its 
solution r is periodic of period T = T{u, c). 

Below, we consider the case c < 0. Then either r < (the left case) or r > (the 
right case). To express r through elliptic functions, we denote by tq > ri > the 
positive solutions of — wp^ — c = 0. Then ri < |r| < tq and one can rewrite (12.30 
as 

r'2 = |(|r| - ri)(ro - \r\){\r\ + Tq + n - fw). (2.4) 

Therefore 2ro + ri > tq + 2ri > Introducing a new variable s G (0, 1) via 

|r| = ri + (ro — ri)s^, we transform (12.41) into 

where a, k, k' are positive constants {k"^ + k'"^ = 1) given by 

^2 _ 4ro + 2ri - 3a; ^2 _ 2ro - 2ri ^,2 _ + 4ri - 3a; 
" ~ 12 ' ~ 4ro + 2ri -3a;' ~ 4ro + 2ri -3a;' 

Therefore 

|r(x)| = ri + (ro — ri)cn^(ax; /c). (2.5) 



a \/uj V A/a; 



As usual, here and below, K{k) and E{k) denote the complete elliptic integrals of the 
first and the second kind in a Legendre form. Let us recall for later use the system 
they satisfy: 

kK' = -K, kE' = E- K. 

1 — k'^ 

Lemma 2.1. For any a; > and T & I , there is a constant c = c{u)) such that the 
periodic traveling-wave solution (12. 5 p determined by H{r,r') = c{u) has a period T. 
The function c{uj) is differentiable. 

Proof. The statement follows from the implicit function theorem. It is easily seen 
that the period T is a strictly increasing function of k: 

'2_T\l^/l.\<')/-,_,2,l.i\r^li 



r.yi — wiTTAKikw - fc(2fc^-i)i^(fc)+2(i-fc2+fc-t)K'(fc) 

dky^ ^ K +KJ\[^K)) - 2(l-fe2+fe4)3/4 



2{l-k^+k'^)E(k)+(l-k'^){k'^-2)K{k) PI 
2fc(l-fc2)(l-fc2+A:4)3/4 > U. 



Given u and c in their range, consider the functions ro(a;, c), ri(cj, c), k{u,c) and 
T{u, c) given by the formulas we derived above. We obtain 



dT dTdk 1 dT d{k 



2^ 



dc dk dc 2k dk dc 
Further, using that fc^ = 3 4y^j^2n -3t^ ^ ^ I^o ~ ~ I'"? ~ '^^i^ have 

dc (4ro + 2ri - Su)^ 

3a;^(ri - ro) 



2' 



We see that dT{u, c)/dc^ 0, therefore the imphcit function theorem yields the result. 
□ 

Equation (12.1 p has the following conservation laws 

Q{u) = i Uxudx, P{u) = / \u\'^dx, E{u) = / {\ux\'^ — )dx. 

Jo Jo Jo 

Let us consider the functional 

M{u) = E{u) + + v^)P{u) - 2vQ{u). 
Next we introduce the pseudometric 

d(u,ip)= inf \\u(x,t)-e'\(x-^,t)\\i. (2.7) 

For a fixed g > 0, we denote 

dliu,cp)= inf (||«,(x,t)-e>.(a;-e,t)|P + g|Kx,t)-eXx-e,t)|p). 

(r,,Oe[0,27r]x[0,T] 

(2.8) 

Clearly, the infimum in (12.81) is attained at some point {r],^) in the square [0, T] x [0, T]. 
Moreover, for q E [qi, ^2] C. (0, 00), (12. 8 p is a pseudometric equivalent to (12. 7p . 

Lemma 2.2. T/ie metric dq{u,ip) is a continuous function oftE [0,oo). 

Proof. The proof of the lemma is similar to the proof of Lemmas 1, 2 in [S] □. 
Now, we can formulate our main result in this section. 

Theorem 2.1. Let (f be given by (IL2p . with r 7^ 0. For each e > there exists 
S > such that ifu{x,t) is a solution of (12. ip and d{u,(f)\t=o < S, then d{u,(f) < e 

Wt e [o,cx)). 

The crucial step in the proof will be to verify the following statement. 

Proposition 2.1. There exist positive constants m,g,5o such that if u is a periodic 
solution of (12. ip . u{x, t) = u{x + T,t), P{u) = P{(p) and dg{u, ip) < 6q, then 

M{u) - M{ip) > md^iu, if). (2.9) 
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Proof. In order to estimate AM = M{u) — M{if), we set 

u{x, t) = e^'^ip{x - ^, t) + h{x, t), 



( = v^x-0 + {u- v^)t + 7], F{s) = -||e> + hs\^ 



We have 



H(-H3 . M3, . f (1, - f (0, ^ no) . ^ . 

where 

< s < 1, F'(0) = -2\ip\Reie'\h), F"(0) = -|^| [Reik'^e-'''^) + 3\h\^] . 
Integrating by parts in the terms containing hx and h^, we obtain 
AM = M{u) - M{ip) 

= 2Re e'''^[-ipxx + {u + v"^ - \<^\)ip + 2ivipx]hdx 

+ + {io + v'- l\ip\)\h\^ - Mi?e(e-2</i2) _ 2tvhK]dx 

= Jl + /2 + /3. 

Using that r{x) satisfies the equation ( 12. 2p we obtain that Ji = 0. 

Let h{x, t) = {hi +i/i2)e*^, where hi and /i2 are real periodic functions with period 
T. Then we have 

|/i|2 = hl + hl 

l^xP = hj, + hi, + 2v{hih2x - hiM + v\hl + hi) 



(2.10) 



i?e(e-2</i2) = hl-h 



Jq hhxdx = i Jq {hixh2 — hih2x — vh\ — vhl)dx. 
Finally, for I2 we obtain 

h = [ [hi, + {co - 2\r\)hl]dx + [ [hj, + {u - \r\)hl]dx = Mi + M2 
Jo Jo 



Consider in [0,T] — [0,2K{k)/a\ the formal differential operators 

supplied with periodic boundary conditions. By the above formulas, ro — ri = Ga'^k'^, 
2ro — uj = 4«^(1 + k"^). Taking y = ax as an independent variable in Li, one obtains 
Li = a^Ai with an operator Ai in [0, 2K{k)] given by 

Ai = + a ^[uj- 2(ri + (ro - ri)cn^{y; k))] 

, - 2ro 2(ro - n) ^ 
- H n sn [y; k) 



4(1 + A;2) + 12A;W(2/;A;). 



dy2 

The spectral properties of the operator Ai in [0, 2i^(/c)] are well known. The first 
three (simple) eigenvalues and corresponding eigenfunctions of Ai are 

/xo = - 2 - 2VT^^WTW < 0, 

jjjoiy) = dn{y; k)[l - (1 + 2k'^ - y/1 - k'^ + '^k^)sn\y; k)] > 



i>i{y) = dn{y; k)sn{y; k)cn{y; k) = \j^cn^{y] k) 
;i2 = A;2 - 2 + 2Vl - F + Ak^ > 

iP^iy) = dn{y; k)[l - (1 + 2k'^ + VT^^k^TW)sn'^{y; k)]. 

Since the eigenvalues of Li and Ai arc related by A„ = a'^fin, it follows that the first 
three eigenvalues of the operator Li, equipped with periodic boundary condition on 
[0, 2K{k)] are simple and Aq < 0, Ai = 0, A2 > 0. The corresponding eigenfunctions 
are ipo{ax) , ipi{ax) = const.r' and ■02(«a;). 

In a similar way, with L2 — q;^A2, one obtains in [0, 2K{k)] 

d"^ 

A2 = --r^ - 2(1 + k^) + 6khn\y; k) + u/2a'. 
dy^ 

To express ui through a and A;, one should take into account the fact that in the cubic 
equation we used to determine ro and ri, the coefficient at p is zero. Therefore, 

roTi + (ro + ri){\uj - ro - ri) = 0. 

As ro = 2a^ + 2a'^k'^ + ^cu, ri = 2a'^ — Aa^k"^ + |a;, after replacing these values in the 
above equation one obtains a;^ = 16q;^(1 — fc^ + fc^). Since a; > 0, we finally obtain 

A2 = --^ + 2(-l - A;' + Vl - A;2 + A;^) + Qk^sn'^iy- k). 
dy^ 
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On the other hand, (12.50 yields 

|r| = 2a^[l + P + Vl - A;2 + A;^ - 3A;W(y; A;)]. 
The first three eigenvalues and corresponding eigenfunctions of A2 are as follows: 

Co = 0, Vq = r, 

ei = 2-e + 2V1 - A;2 + k^, vi = dn'{y; k) 

€2 = 4 VI - F + k\ V2 = l + k^- Vl-k'^ + - Sk^sn\y; k). 



Estimates for M2. 

From the above explanations we know that when considered in [0, T], the operator L2 
has an eigenfunction r corresponding to zero eigenvalue and the rest of the spectrum 
is contained in (a^Ai, 00). 

In the formulas which follow, we take r = r{x) with an argument x = x — ^ — 2vt. 
The values of ^ and 77 are chosen so that the infimum in (12.81) is attained at that point. 
Therefore the derivative of dg{u, ip) with respect to t] is equal to zero. Together with 
(O, this yields 



I h^e ^^(fx — (fxC^^hx + q{he "^^(f — ^pe^'^h)dx 



2Im {-(fxx + qip)e^'^hdx 



(2.11) 



lo [(^^ — 00 + q + \r\) r/i2 + 2vr'hi]d. 



X. 



We set h2 = /3r{x) + 6, 6rdx = 0. Substituting in (12. lip , we obtain 



Using that 



|r3/2||2 



,3/2112 



V — UJ + q + 



r 

I [6r\r\ + 2vr'hi\dx = 0. 
Jo 



> UJ (see estimate Al of Appendix 1), we obtain the estimate 



< 





r 


+ 2vr'hi)dx 




\\r^\\ ■ \\d 


+ g)||r| 
1 +2|t^| 




\\hi\\ 



{v"^ + g)||r| 
< modioli + ||/ii||). 
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where ttlq = 2mi {v,u)/{q + u"^) and 



^21 



\v\\ 




\r 


\r ■ 


- U)r\\ 









mi(f , Lo) = max 

ce[-§a;3,0] 

Clearly, the first and the third terms are uniformly bounded for u fixed. The bound- 
edness of the second and the fourth ones follows from the estimates in Dl of Appendix 
l.(The third and fourth terms are included for later use.) We will use below that for 
V and u fixed, mo — when g — )■ oo. Further, 

\\h2\\ < |/3|||r|| + ||e|| < modioli + ||/ii||) + 11^^11 = (mo + 1)||^|| + mo||/ii||. 



Hence, we obtain 



2(mo + l)2 V"^o + l 



\\h\\\ (2.12) 



Since = and {6,r) = 0, then from the spectral properties of the operator L2, it 
follows 

M2 = {L2h2,h2) = {120,9} > ahMl''. 
From here and (12.121) . one obtains 

2(mo + 1)^ (mo + I)'' 

Estimates for Mi. 

We set 

^1 = 7i^o(^) +72r'(x) + 6'i, r{x) = utpo{x) + (2.14) 

where 

(^1, ^0) = (^1, r') = V-o) = (V-o, r') = r') = (2.15) 
and 7i, 72 and z/ are some constants. By ( I2.15p . we have 

M,{h,) = {L,h,,h,) = 7i'Ao(^o,^o) + {Lie,,e,). 

Therefore, from spectral properties of the operator Li it follows 

Mi(/ii) >7?Ao||^o||' + A2||^i|p. (2.16) 

The fundamental difficulty in the estimate of Mi is the appearance of the negative 
term 7^Ao| iV'ol Below, we are going to estimate it. From the condition 

P{u)= f \h + e'^'^^{x-i,t)\^dx = P{^) 
Jo 

we obtain ^ ^ 

\\h\\^ = 2Ref e''^''^(x - ^,t)hdx = -2 I rh^dx. 
Jo Jo 
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Then using fl2.14p . we have 

-^l|/i||' = ^7ill^oll' + ^ i^e^dx 

and therefore ^ 

iim? = ^T^Alm? + r ^e,dx\ . (2.17) 







From (I2.17p . we obtain 



iim? < (ii^iiAir + '!^urm?) , (2.18) 

where c? is a positive constant which will be fixed later. Using fl2.17p and fl2.16p . we 
derive the inequality 



Below, we will denote by C^,, Dm positive constants. By using (12.181) and (I2.19p . 
we derive the inequality 



>CiA2ii^ir-i^ 



(see the estimates in point CI of the Appendix 1). 

We denote ^ = hi - 72r'(x) = 7iV'o(^) + ^i- Then from (121^]) . (E^O]), we have 



<C2||0l|P + D2|"^"' 

Then 



and hence, by ([2201), 



C2 C2|a|2 



M,>^||tf|p-^i^2%t^||ft||^ (2.21) 



After differentiating ( 12. Sp with respect to ^ and using ( 12. 2p . we obtain 

= 2i?e e'-'^''{(pxxK + q^Pxh)dx 

= 2Re Jq [{r" + 2ivr' — v'^r){hix — ih2x — ivhi — f /12) + q{r' + ivr)(hi — ih2)]dx 
= 2 [{-00 + 2\r\ + 3f ^ + q)r'hi + v{v'^ - 3oo + 3\r\ + q)rh2]dx. 
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From (12. lip , we have 



qrh2dx = — / [{v — u + \r\)rh2 + 2vr'hi]dx 
Jo 



and replacing in the above equahty, we obtain 



[{-u + 2\r\ + q)r'hi + v{2\r\ - 2u)rh2]dx = 0. 







Substituting hi = 72r'(x) + in the above equahty and using the orthogonahty 
condition (r', i?) = (r', •jiipo + 6*1) = 0, we obtain 



M /m2 I 2 2||v/fr'||2 

72||r||M-^ + ^' + g+ "m.|2 



[ [v{2\r\-2uj)rh2 + 2\r\r'i^]dx = 0. 
Jo 



As il^ip^^ > CO (see estimate Bl from Appendix 1), 
l72|||r'|| < 





-2uj)rh2 + 2\r 


r''d]dx 


{v^ + q)\\r'\ 





< 2 



\v 




\r 


r — ur\ 


\-\\h2\ 


\ + 2\\\r 


r'\ 


\-\m 


(f^ + q)\\r'\ 





<m,m\ + \\h2\\). 

Hence 

||/^i|| < I72I llr'll + Pll < (mo + 1)P|| + mo||/i: 

which yields 



2 



2 / _ \ 2 

2| 



2(mo + l)2 V"^o + l 
Replacing in (12.2ip . we finally obtain 

>- C^2^ll'-ll= - C^^W.r - '-^^^^^M'. (2.22, 

The estimate for AM. 

From fl2.13p and fl2.22p . we have (fixing q large and therefore tjiq small enough) 

A# _L A/f > CiA2-2C2a^':im§ \ m, ||2 , / C2a^<:i-2Ci Aamg \ ,,, 
_ C.A.D^+C2D. ||^||4>^^||;^||2_^^||,||4 



|2 
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On the other hand, estimating directly I2 from below (for this purpose we use its 
initial formula), we obtain 

h >\\h.\? + lo {u + v^- 2\y,\)\h\^dx - 2\v\ \h\ ■ \hMx 
> WKW^ + {io + v'' - 2max\r\)\\h\\^ - 2v^\\h\\^ - IWKW^ 



= l\\hx\\^ + (cu-v'-2ro 
Let < m < |. We have 
AM = 2ml2 + (1 - 2m) (Ml + M2) + I3 



> m\\h,\\^ + 2m{u -v^- 2ro)||/i|p + (1 - 2m)(C3| P - D^WhW^) + h 



> m\\K\\' + [{2m{uj -v^- 2ro) + (1 - 2m)Cs]\\h\\^ - (1 - 2m)D3||'i|| - M3|- 
We choose m, so that 2mq = (1 — 2m)C-i + 2m{uj — f ^ — 2ro), i.e. 

C3 

g + C3 - w + t.2 + 2ro ■ 
From the continuity of \z\ and |z|e~^*"^'^^, we have 

From the inequality 

\h\^<^j^ \h\^dx + 2(^j^ \h\^dx \K\^d:^' 

we obtain ^ ^ ^ 

\h\'<^[ \h\'dx + ^[ \h\'dx + ^[ IKl'dx. 
T Jo Jo Vq Jo 

Hence for sufficiently large q, we obtain 

2 

max|/i(x, t)p < —-d'^{u,ip) 

and moreover \\h\\^ < q~^dg{u,ip). Consequently we can choose 60 > 0, such that for 
dq{u,ip) < (5o, we will have [max(4|a|5 + |/ip) + (1 — 2m)L'3|a|^]||/i|p < qm. 

Finally, we obtain that if dq{u, ip) < 60, then AM > mdq{u, (p). Proposition 12. II is 
completely proved. □ 



Proof of Theorem 12.11 We split the proof of our main result into two steps. We 
begin with the special case P{u) = P{^p)- Assume that m,q,6o have been selected 
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according to Proposition 12.11 Since AM does not depend on t, t G [0, oo), there exists 
a constant / such that AM < l(P{u,(p)\t=o. Below, we shall assume without loss of 
generality that I > l,q > I- 
Let 

e > 0, = mm 
and d[u,(f)\t=Q < 5. Then 



IqJ 2 ' V / 



dqiu,(p) < q^/'^d{u,(p)\t=o < y 

and Lemma [2l2] yields that there exists a to > such that dq{u,ip) < Sq ii t ^ [0,to)- 
Then, by virtue of Proposition 12.11 we have 

AM > mdliu.if), t G [0,to)- 

Let tmax be the largest value such that 

AM >mc/J t G [0,t„a,). 

We assume that tmax < oo. Then, for t G [0,tmaa;] we have 

, AM / / . 51 

dliu^if) < < -d\u,ip) t=o < -5' < -J- 

^ mm m 4 



Applying once again Lemma 12. 2t we obtain that there exists ti > tmax such that 

dg{u,(p) <5o, t e [0,ti). 

By virtue of the proposition, this contradicts the assumption tmax < oo. Conse- 
quently, tmax = OO, 

AM >mdl{u,ip) >md'^{u,ip), te[0,oo). 

Therefore, 

d (u, if) < < —S^ <e^, te 0, oo), 

m m 

which proves the theorem in the special case. 

Now we proceed to release the restriction P{u) = UmIP = = P{f)- We have 



(16aVl - P + fc4 \{k^ - 2 + Vl - A;2 + k^)K{k) + 3E{k) 



,1/2 



Below, we are going to apply a perturbation argument, freezing for a while the period 
T and the parameters u,c in (12. 3p . We claim there are respective parameter values 
uj*,c*, and corresponding (p*, r*, a*, k*, see (12. 2p . (12. 3 p and (12. 5p . such that ip* has a 
period T in a; and moreover, ||v?*|| = ||^^||- By (12. 9p . we obtain the equations 

2K(k*) ^ 
^ - T = 0, 

a* (2.23) 
||r*|P - = 0. 
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Moreover, one has \ \(p*\ \ = \ \u\ \ and we could use the restricted result we established 
above. As k = k*(T, \ a = a*(T, \ it remains to apply the imphcit function 
theorem to (12.231) . Since the corresponding Jacobian determinant reads 



d 
dk* 



a 

dk* 



2K{k*) 



-* 1 12 



a ( 2K{k*) 

aa* \ a* 



a 

Qa* 



112 



> 



the needed properties are established. 



50000- 



20000 
10000 



Figure 1: Graph of the Jacobian determinant 

By (12.231) and our assumption, we have 

K{k) _ K{k*) _ T 
a a* 2 

Next, choosing r] = 2{uj* — uj)t, ,^ = 0, we obtain inequality 

d\ip\ip) < (l + t;')||r* -r|p + ||r*' -r'ip. 



(2.24) 



Denote for while <I>(p) = 2a^(l - 2A;^ + Vl - fc^ + k^ + Sk'^cn^{ax; k)), where k = k{p) 
is determined from K{k) = \pT. Then r* — r = $(«*) — $(a) = {a* — a)$'(p) with 
some appropriate p. Moreover, |$'(p)| < Cq with constant independent of the values 
with * accent. Hence \r* — r\ < Cola* — a\. Similarly \r'* — r'\ < Ci|a* — a|. All this, 
together with 02.241) yields 

Of 2C 
d{v\ ^) < C\a* - a\C = —\K{k*) - K{k)\ = —\K'{k)\\k* - k\. (2.25) 



Let e > 0. From the inequalities 
it follows that 

-\\r\\6<{M 



Iv^lll < d{u,(f)\t=o < S 
\lP*\\ - 1 < ||r||5 
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and, consequently, | 1 1?"*! P — | l^^l P | < | kl pf^ij where (5i = (1 + 1 1^1 — 1. 
On the other hand, we have (using f l2.24p again) 



1^*1 1 2 1 1 ^ 1 1 2 



2m_ 

7-3/2 



-{K{k)f/^^l-k^ + k^ [{k^ - 2 + VI - A;2 + k^K{k) + 3E{k)\ \ 
> C2\\k*-kl 

(2.26) 

Thus combining fl2.25p and (12.261) . we get 

d{u, ip*)\t=o < d{u, if) \t=Q + rf(v9, if*) \t=Q <5+\ |r I pC^i = 5o- 

We select 5 sufficiently small and apply the part of the theorem which has been 
already proved, 

ci(M,v2*)|t=o < ^0 ^ d{u,ip*) < |, t G [0, oo). 
Choosing an appropriate 5 > 0, we obtain that 

d{u, if ) < d{u, if*) + d{ip, if*) < e, 
for all t G [0, oo). Theorem 12.11 is completely proved. □ 

3 Existence and stability of periodic traveling waves 
for the cubic Schrodinger equation 

Consider the cubic nonlinear Schrodinger equation 

iut + Uxx + \u\^u = (3.1) 

where u = u{x,t) is a complex-valued function of (x,t) G M^. 

Equation (13.11) possesses the following family of traveling-wave solutions 

if{x, t) = e^(''^+('^-^')*)r(x - 2vt), (3.2) 

where u and v are real parameters and the real- valued function r{x) satisfies equation 

r" -ur + r^ = 0. (3.3) 

Integrating once again, we obtain 

r'2 - oor^ + lr^ = c (3.4) 

14 



and r is a periodic function provided that the energy level set H{x,y) = c of the 
Hamiltonian system dH = 0, 



Hix,y) = -ux^ + ^a;^ 

contains an oval (a simple closed real curve free of critical points). The level set 
H{x,y) = c contains two periodic trajectories if w > 0, c G (— |a;^,0) and a unique 
periodic trajectory if w G M, c > 0. Under these conditions, the solution of fl3.3p is 
determined by H{r, r') = c and r is periodic of period T = T{u, c). 

Below, we are going to consider the case c < 0. Let us denote by ro > ri > 
the positive roots of — ur"^ — c = 0. Then, up to a translation, we obtain the 
respective explicit formulas 

r[z) = ^rodniaz; k), k = ^ — = 5 -, a = —=, T = . (3.5) 

Recall that K{k) and E{k) are, as usual, the complete elliptic integrals of the first 
and second kind in a Legendre form. By (13.51) . one also obtains u = {2 — k'^)a^ and, 
finally, 

T^'-^™^. ^^(0.1), T.,.{^^^. (3.) 



We take w > 0, an energy level c G (— |a;^,0) and let T be the (minimal) period 
of r{x). Assume that f 7^ is chosen to satisfy vT/2Ti G Z. Then there are two 
periodic solutions r of (13.31) : left (r < 0) and right (r > 0) and the corresponding 
traveling-wave solution Lp{x,t) is periodic with respect to x of period T. 



Lemma 3.1. For any > and T ^ I, there is a constant c = c{u)) such that the 
periodic traveling-wave solution (13. 5p determined by H{r,r') = c{u) has a period T. 
The function c{uj) is differentiable. 

Proof. See [TO], Lemma 3.1. 

Equation (13.11) has the following conservation laws 

Q{u) = i Uxudx, Piu) = / \u\'^dx, E{u) = / {\ux\'^ 7^)dx. 

Jo Jo Jo 2 

Let us consider the functional 

M{u) = E{u) + (w + v'^)P{u) - 2vQ{u). 

Theorem 3.1. Let if he given by (13. 2p . with r ^ 0. For each e > there exists 
6 > such that ifu{x,t) is a solution of (13. ip and d{u,(f)\t=o < then d{u,(p) < e 
\/t G [0,00). 

The crucial step in the proof will be to verify the following statement. 
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Proposition 3.1. There exist positive constants m,q,SQ such that if u is a periodic 
solution of (13. ip . u{x, t) = u{x + T,t), P{u) = P{(p) and dg{u, ip) < 60, then 

M{u) - M{ip) > md]{u, (^). (3.7) 

Proof. In order to estimate AM = M{u) — M{if ), we set 

u{x, t) = e^^'ip{x — ^,t) + h{x, t), 

C = v{x-0 + {^ + v^)t + r], 
and integrating by parts in the terms containing h^ and h^, we obtain 

AM = M{u) - M{ip) 

= 2Re Jq e'''^[—(pxx + (w + f ^ — Iv^Hv^ + 2iv(px]hdx 

+ lo [l^^l^ + {u} + v^- 2|^|2)|/i|2 - |(^|2i?e(e-2»C/i2) _ 2ivhK]dx 

-\ |/i|2(4i?e(e^V^) + \h\^)dx 



= h + l2 + h. 

Using that r{x) satisfies the equation (13. 3p we obtain that Ji = 0. 
Let 

h{x,t) = {hi + ih2)e'^, 

where hi and /12 are real periodic functions with period T. Using (I2.10p . for I2 we 
obtain 



[hl^ + {u - 3r^)hl]dx + / [hl^ + {u - r^)hl]dx = Ml + M2 

Jo 

Introduce in L'^[0, T] the self-adjoint operators Li and L2 generated by the differential 
expressions 

Li = -dl + {uj-3r^), 

(3.8) 

^2 



Lo = -dl + (00- r') 



with periodic boundary conditions in [0,T]. 

We use now (13. 5 p and (13.61) to rewrite operators Li, L2 in more appropriate form. 
From the expression for r{x) from (13. 5 p and the relations between elliptic functions 

sn{x), cn{x) and dn{x), we obtain 

Li = a^[-dl + Qk'^sn^{y) - A - k"^] 
16 



where y = ax. 

It is well-known that the first five eigenvalues of Ai = —dy + 6k^sn^{y,k), with 
periodic boundary conditions on [0,4:K{k)], where K{k) is the complete elliptic inte- 
gral of the first kind, are simple. These eigenvalues and corresponding eigenfunctions 
are: 



^0 


= 2 + 2P -2^1 


- F + k\ 


My) = 


1 - (1 + P - Vi - 


-k^ + k^)sn\y,k), 


1^1 


= l + k^ 




My) = 


cn{y, k)dn{y, k) = 


sn'{y,k), 




= 1 + 4P, 




My) = 


sn{y, k)dn{y, k) = 


-cn'{y,k), 




= 4 + P, 




My) = 


sn{y, k)cn{y, k) = 


—k^'^dn'{y, k), 


U4 


= 2 + 2A;2 + 2Vl 


-k^ + k\ 


My) = 


1-ii + k^ + Vi- 


-k'^ + k^)sn^{y,k). 



It follows that the first three eigenvalues of the operator Li, equipped with periodic 
boundary condition on [0,2i^(A;)] (that is, in the case of left and right family), are 
simple and Aq = a^(z^o — ^'s) < 0, Ai = a^(z/3 — ^3) = 0, A2 = 0^(1/4 — z/3) > 0. The 
corresponding eigenfunctions are ipQ = (f)o{ax),^pl = r'{x),'ip2 = (^^{ax). 

Similarly, for the operator L2 we have 

L2 = + 2k'^sn^{y, k) - k'^] 

in the case of left and right family. The spectrum of A2 = —dy + 2k'^sn'^{y, k) is formed 
by bands [fc^, 1] U [1 + A;^,+oo). The first three eigenvalues and the corresponding 
eigenfunctions with periodic boundary conditions on [0, 4K{k)] are simple and 

eo = /c^ 9o{y) = dn{y,k), 
ei = 1, Oiiy) = cn{y,k), 

e2 = l + k'^, e2{y) = sn{y,k). 

From (13.41) it follows that zero is an eigenvalue of L2 and it is the first eigenvalue 
in the case of left and right family, with corresponding eigenfunction r{x). 

Estimates for M2. 

As in Section 2, we use below r = r{x) with an argument x = x — C, — ^.vt. From 
the explanations above, we know that when considered in [0,T], the operator L2 has 
an eigenfunction r corresponding to zero eigenvalue and the rest of the spectrum is 
contained in (a^, 00). 

The derivative of dg{u,(p) with respect to 77 at the point where the minimum is 
attained is equal to zero. Together with (13. 3p . this yields 

T 

[{v"^ - u + q + r'^) rh2 + 2vr'hi] dx = (3.9) 
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We set h2 = /3r{x) + 6, 6rdx = 0. Substituting in (I3.9p . we obtain 



/3||r 



|2 ..2 



f — w + g + 



U2||2 



r 



12 



+ /" [9r^ + 2vr'hi]dx = 0. 
Jo 



Using that > w (see estimate A2 of the Appendix2), we obtain the estimate 

jQiOr^ + 2vr'hi)dx 



r < 



(g + f2)||r 

I I T"^ I I 

< 



1 1 1 


|-||^ll + 2|t; 


r 


'\\-\\hi\\ 




(g + f^) 


|r| 





<mom\ + \\hi\\).. 
where mg = 2mi (t>, (X')/(g + w^) and 

,,^311 oL,l ll^'ll oL,l 11^3 , ,^\\ Qll»,2^/| 

mi(f , uj) = max 



ce[-iaj2,0] 



5 



2 






^3 _ 











(the third and fourth item are included for later use). It is obvious that the first and 
the last fractions are bounded. For the second and the third ones, see estimates D2 
in Appendix 2. We will use below that for v and u fixed, rriQ — )■ when q — t- oo. 
Further, 

\\h2\\ < m \\r\\ + \\e\\ < modioli + WhW) + \\e\\ = (mo + 1)||^|| + mo||/ii||. 
Hence, we obtain 

m\' > - ( ii/^iir. (3.10) 

" " - 2(mo + l)2 V^o + i; 
Since = and {6,r) = 0, then from the spectral properties of the operator L2, it 



follows 



M2 = {L2h2,h2) = {126,6) > a^6,6) > 



From here and fl3.10p . one obtains 



2 



Estimates for Mi. 

We set 

hi = 'ji^o{x)+'j2r'{x) + 6i, r{x) = utpo{x) + (3.12) 

where 

(01, ^0) = {61, r') = (tA, V-o) = (^0, r') = (^, r') = (3.13) 
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and 7i, 72 and v are some constants. As the same line as for Mi, we obtain 

1 



T X 2 



7?ll^oir = ^^(k/^lP+ / ^Mx) . (3.14) 



and 

^1 >(A2 + Ao(l + i)45|^)||^i|P + ^||/.||^ 

(3.15) 

(see the estimates in point C2 of the Appendix2). 

We denote 'd = hi — 72r'(x) = 71-^0(2;) + ^i- Then from fl3.14p . f l3.15p and (since 
A2 = —OJ + V 4a;^ + 6c and A2 = — — A/4a;^ + 6c the inequahties A2 < ^|Ao| < w, we 
have 



4,.2||^0||2 



Then 

ii»,ip>"''"^ ft"-' 



C2 (720; 2 



and hence, by f l3.15p and A2 < w. 



^1-^2:10112 C1D2 -\- C2D1 i||,||4 



O2 L/2 
= C3A2||^|P-/^3C^5||/,||4. 

After differentiating (12.80 with respect to ^ and using (13. 3p . we obtain 

= 2Re e''^{(pxxhx + q(pxh)dx 

= 2 [{-Lo + 3f ^ + 3r2 + g)r'/ii + v{-3uj + v'^ + 3r^ + q)rh2]dx. 
From (13.90 . we have 



(3.16) 







qrh2dx = — [2vr'hi + (— a; + + r'^)rh2]dx 
Jo 

and replacing in the above equality, we obtain 



f [{-u + + 3r2 + qyhi + v{-2uj + 2r'^)rh2]dx = 0. 
Jo 

Substituting hi = 72?"' (x) + {} in the above equality and using the orthogonality 
condition {r',{}) = {r',jiipo + 61) = 0, we obtain 



72||r'||^ [-uj + v^ + q + 



3 rr' 


2 






2 



^ + ^ [2f (-w + r^)r/i2 + 3rV^9]dx = 0. 
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Using that ^ > Lo (see the estimates in point B2 of the Appendix2), we further 
have 



Hence 

which yields 



I72|||r1| < 



X 



{v^ + q)\\r'\ 



< 



2 


V 




— ur + r^\ 


\-\\h2\ 


1 +3||rV| 


\-\m 






{v^ + q)\\ 







<mo(P|| + 11/^211). 
||/^i||<|72|||r'|| + P||<(mo + 



+ mo\\h2 



2> ll^ilP 



2(mo + l)2 V"^o + l 
Replacing in fl3.16p . we finally obtain 

,2 



Ml > 



^3-^2 



2(mo + ly 



l|, ||2 C3X2mQ 



(3.17) 



The estimate for AM. 

From f l3.1ip and fl3.17p . one obtains 



Ml + Ma > 



2(mo + l)^ 



4(mo + 1)^ 



We now fix q so that uml < IC3A2 and assuming that C3 < | (which is no loss of 
generality), one has also 4C3A2mo < \uj. Therefore we come to 



M1 + M2 > CMWhW^ + \\h2\\^) - Dstu^^ 



C4X2 



where C4 and D3 are absolute positive constants independent on the parameters of 
the system. 

On the other hand, estimating directly I2 from below (for this purpose we use its 
initial formula), we have 



h > + {u + - 2r^)\h\^dx ~ 2\v\ \h\ ■ \h,\dx - r^\h\^dx 



> + {io + ~ 2rl 



2v'\\h\\'-m^ 



l\\h,\\^-{v' + 5co)\\h\\\ 
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Similarly, I/3I < max(4|2w| 2 + Let < m < |. We obtain 

AM = 2ml2 + {1 - 2m){Mi + M2) + I3 

> m||/i^||2 - 2m(t;2 + 5uj)\\h\\'^ + (1 - 2m){CA2\M'^ - ^s^^^ll^l 
-max{A\2u\^h\ + \h\'^)\\h\\'^ 
= m||/i^||2 + [-2m{v^ + 5uj) + (1 - 2m)Ci\2] 
-[max{4\2uj\^h\ + j/ip) + (1 - 2m)D3UJ^\h\\'^]\\h 



2 



We choose m, so that 



C A 

2qm = (1 - 2m)CiX2 - 2m(v^ + 5u), i.e. 2m = ^ — — — — — < 1. 

From the inequality 

T T T ^ 

\h\^<^ \h\^dx + 2( \h\^dx \h^\^dx 
T Jo \Jo Jo 

we obtain ^ ^ ^ 

\h\^ < I \h\^dx + y/ql \h\^dx + —l \h^\^dx. 
T Jo Jo \/Q Jo 

Hence for sufficiently large q, we obtain 

2 

max|/i(x, t)P < dl(u,ip) 

and moreover \\h\\^ < q'^d'^{u,ip). Consequently we can choose 60 > 0, such that for 
dq{u,(p) < 5o, we will have [max(4|2w|5 + \h\'^) + (1 — 2m)D3U^\\h\\'^ < qm. 

Finally, we obtain that if dq{u, ip) < 60, then AM > mdq{u, ip). Proposition 13. II is 
completely proved. □ 

Proof of Theorem 13.11 The proof of the theorem in the case P{u) = = 
llv^lP = Pi^p) is the same as in Theorem 12.11 If P{u) 7^ Pi'p), we proceed similarly 
as in Theorem 12.11 We have \\(p\ \ = {2\/2roE(k)y^'^ , where ro is given by (13. 5p . 

We claim there are respective parameter values u}*,c*, and corresponding (p*,r*,rQ, 
k*, see (13. 5p and (13. 6p . such that (p* has a period T in x and moreover, 2\/2rQE{k*) = 
By (13. Sp . we obtain the equations 

2V2Kik*) ^ _ ^ 

ro ' (3.18) 

2V2r*E{k*) - = 0. 

If (I3.18P has a solution k* = A;*(T, | |m| |), Vq = rQ(T,\\u\\), then the parameter values 
we need are given by 

2w* = (2 - k*yf, 2c* = {k*^ - l)r*^ 
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Moreover, one has \ \(p*\ \ = \ \u\ \ and we could use the restricted result we established 
above. As k = k*(T, \ \'^\\), tq = r^iT, \ \'^\\), it remains to apply the implicit function 
theorem to (13.181) . Since the corresponding Jacobian determinant reads 



2V2K'(k*) 2y/2K(k*) 
r* V*^ 

2V2r*E'{k*) 2^E{k*) 



-{KE)' > 



the existence of u* and c* with the needed properties is established. 
By (13.181) and our assumption, we have 

2V2Kik*) ^ 2V2K{k) ^ ^ 

Choosing r] = 2{uj* — uj)t, ,^ = 0, we obtain the inequality 

(f{ip\^) < {l + v^)\\r* -r\\^ + \\r*' -r'W^ 

For = pdn{zp, k{p)), where k = k{p) is determined from K{k) = \pT, we 

get r* — r = ^{tq) — $(ro) = {tq — ro)$'(p) with some appropriate p. Moreover, 
W{p)\ < Co. Hence \r* - r\ < Co|r* - ro|, \r*' - r'\ < Ci|r* - ro|, and by ( l3l9l) 

d{v\ < \r* - rol < ^^\Kik*) - Kik)\ < '^-^\K\k)\\k* - k\. 
From the inequalities 

lllv^ll - llv^lll = 111^*11 - llv^lll < d{u,(f)\t=Q < 6 

it follows that 

- (2V2roEy^^^ 6 < (| |^| |)"'| ly^l | - 1 < (2V2roE'^'^^^ 



6 



and, therefore, 1 - 6i < """^^^^ < 1 + 5i, i.e. \rQE{K*) - vqEI < tqESi, where 
(5i = (1 + (2ro^)-^/25)2 - 1. On the other hand 

\r;E{k*)-roE{k)\ = '^\K{k*)E{k*)-K{k)E{k)\ = '^\{KEy{k)\\k*-k\ > C2\k*-k\ 

with appropriate C2 > 0. Thus 

d{u, (p*)\t=o < d{u, ip) \t=o + d{ip, ip*)\t=o <S + roEC6i = 60- 
The rest of the proof is the same as in Theorem 12.11 □ 
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4 Appendix 1 



Below, we provide some estimates needed in our proofs concerning the quadratic 
Schrodinger equation. Without loss of generality, we will assume that r > 0. The 
case of negative r is dealt with by changing r — — r in all equations from Section 1. 

For n e Z and c G {—^u^, 0), consider the line integrals /n(c) and their derivatives 
4(c) given by 

/„(c) = / x-ydx, I'M = / ^ (4.1) 

Jh=c J h=c ■^y 

where the integration is along the right oval contained in the level set {H = c} and 
H{x,y) = y"^ — ujx^ + These integrals would be useful because 

r^{t)dt = 2 / r"(t)dt = 2 / = i ^ = 21'^. (4.2) 

Jo Jri Jc + UJx'^-lx^ Jh=c y 

(we applied a change of the variable r(t) = x in the integral and used equation (12. 3p ). 
The properties of /„ are well known, see e.g. [9] for a similar treatment. Below, we 
list some facts we are going to use. 

Lemma, (i) The following identity holds: 

3ncln-i + 3u{n + 3)/„+i - {2n + 9)In+2 = 0, n eZ 

which implies 

h = coh, h = ^c/o + f^uj-'h, h = ^^culo + (^c + iiu;=^)/i, (4.3) 

(ii) The integrals Iq and Ii satisfy the system 

GcTq + 2^2// = 5/^^ 

Gojcl'o + (30c+ 12cj3)j/ ^ 35j^_ 

(iii) The ratio R{c) = I[{c)/rQ{c) satisfies the Riccati equation and related system 

c = 6c(3c + w^), 
R = oj^R^ + 6cR - 3uc 



c = Qc(3c + u'^) 

6c(3c + uj^)R'{c) = oo^R^ic) + 6cR{c) - 3wc, . ^ (4.4) 



which imply estimates 

- — < Ric) < — -—. 4.5 
6 2uj^ 

The equations in (i)-(iii) are derived in a standard way, see [9] for more details. 
The estimates (14. 5 p follow from the fact that, in the (c, i?)-plane, the graph of R{c) 
coincides with the concave separatrix trajectory of the system (14. 4 p contained in the 
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triangle with vertices (0,0), (— |cj^,cj) and (0,^) and connecting the first two of 
them. 

After this preparation, we turn to prove the estimates we used in the preceding 
Section 2. 

Al. The estimate for Ai = ^-^j^- By (g^]), dMD and the first inequality in 
(14. 5p . we have 

^ _ lo ^^^^ _ 4 _ 3 c/^ + 2^2// ^ 3c 6u; ^ ^ 



Bl. The estimate for Bi = By (Q and (jOD, we have as above 



||r'||2 

2 J^r{c + ur^ - y)dt _ 2{cl[ + ul'^ - 



(c + u;r2 - lr^)dt cl'^ + ujI'^ - f/^ 



3 

_ 6 cw/^ + (5c + 2u^)I[ _ 6 cw + (5c + 2uj^)R ^ 12 
~ 7 3c/^ + ~ 7 3c + " Y'^' 

To obtain the last inequality, we used both estimates in (14. 5p . 

CI. The estimate for Ci = A2 + Ao(l + ^) J'ffp - By (1^]) and we have 



|r| I — I |'?/'o| I , where v 



llV^o 



12 ■ 



Therefore 

— ^-^)(^- 

Below, we need to use the following well-known equalities (see e.g. [7j) which are 
written out here for reader's convenience: 

i) dn^xdx = E{k) 

ii) /(f dn^xdx = ^E{k) + ^K{k) 

iii) f^"" dn^xdx = 8fc^-23fe^+23 ^(fe) + 4(/c2-lK2~fc^) ^(fe) 

iv) /o"^ dn^xsn^xc/x = ^E{k) + ^K{k) 

iv) rfn^xsnSrfx = 8^1^^E(A;) + H(i^^!ili±H^K(A:) 

v) J^"" cn^xdx = ^[E{k) + (F - l)K{k)] 

vi) sn^xc/x = ^[i^(A;) - E{k)] 



24 



(vii) dnxdx = J^^ dnxsn^xdx = J^^ dnxcn'^xdx = it /2 

(viii) J^^ dnxsn'^xcn'^xdx = tt/S 
J^^ dn^xdx = ^-^TT 



Now, we first recall that 

r{x) = a^[l + A;2 + - + - 3khn\ax; k)] 



M^) = dn{x; k)[l - (1 + 2P - - k^ + Ak^)sn^{x; k)] > 0. 
From (12.21) . we obtain after integration 

||r|p = w / r{x)dx = wriT + w{rQ — ri) / cn'^{ax)dx. 
Jo Jo 

Then, by using (v) and the expressions of w, ri etc, we calculate 



||r||2 = 16aVl- A;2 + fc4p2 _ 2 + Vl - P + k'^)K{k) + 3E{k)] 
Similarly, by direct calculations and making use of (i)-(ix), we come to the expressions 

A ./] _ t.2 _|_ AM 

llV^oll^ = -1 ±^ {(A;2 _ i)[(8A;4 + 3k^ + 2) - 2(2^ + 1) Vl - F + 4A;4]ir(A;) 



+ [{8k^ - 3P - 2)Vl -k^ + Ak^ - 2{8k^ - 4k^ - P - l)]E{k)}. 



{r,ipo) = vra[(l + k"^ + Vl - k^ + k^){l - 2k'^ + Vl - F + 4^^) 
- — (- - 2A;2 + Vl-P + 4A;4)] 
Next, calculating the asymptotical expansions near = 0, we obtain 



r 



|2 = 167ra3(l-|P + |A;4 + ...) 



(r,^o) =47ra(l-|F + gA;4 + ...) 



llV^olP = -Vl-P + 4A;4(1 - + i|A;4 + . . .) 



and, finally. 



A2 - Ao VI - A;2 + 4fc4 



A2-A0 (r,V^o)' 32 
Let us denote by 1 — 5 the right-hand side of (14.61) . Clearly, 5 = 5{k) satisfies 1 > 5 > 
for k small enough, say < k < ko. Fixing such a small k, then (14. 6 p yields 
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Therefore 



Ci = A2 + Ao 1 + 



d 



- 1 1 > A2 



l+(l + ^)(5-l) 



CA2 



where C is a positive constant, provided that d is chosen sufficiently large. Note 
that the above estimate is not uniform in k when k tends to zero. This is because 
S = 0{k^) and hence d = 0{k-^). 

When k > ko, we can simply draw the graph of the corresponding function f{k) = 



1 + ^ to see that it is negative and placed far from zero. Hence, Ci > CA2, 



too (uniformly for k > ko). 



Figure 2: Graph of f{k) < 



Dl. The estimates for D11 = ^^^m .m''^^ and £'12 = ttit- Making use of statements 

ii llr'll 1^ ||r|| 6 

(i) and (ii) of the Lemma, we have 



_ 5 -3cuI'o + (6c + u^)I[ _ -ujIo + 2/i 

~7" Scl'o + ujH[ ~ To ~ ^' 

The last inequality follows from the fact that the periodic trajectories of the Hamil- 
tonian system lie inside the saddle loop, in the domain H{x,y) < which implies 

< x^{uj — |x) and x < 

Similarly, 

J^{c + u;r'-lr')dt _ci;, + ujl!,-ll'. 



- 



'12 



I'o 



2 
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The last inequality is due to the known fact [H] (Theorem 12), that the function 
Q{c) = /i(c)//o(c) is strictly decreasing, therefore Q(c) > Q{0) = ^u, the value 
obtained by direct calculation of elementary integrals. 



5 Appendix 2 

First of all, let us mention that the operator Li, defined by differential expression 
(13. 8 p and equipped with periodic boundary conditions in [0,T], has the following 
spectral data 

Ao = —uj — V4:u'^ + 6c, ^0 = 3r^ — 30; — Aq, 

Ai = 0, ^i = r', (5.1) 



A2 = —CO + v4w^ + 6c, 1^2 = 3r^ — Sou — A2. 

For n G Z and c G (— |w^,0), consider the line integrals /„(c) and their derivatives 
j;(c) given by 

f f dx 

4(c) = i x''ydx, In{c) = i — — (5.2) 

J H=c J H=c 

where one can assume for definiteness that the integration is along the right oval 
contained in the level set {H = c}, H{x,y) = y'^ — u;x^ + \x^. As above, these 

integrals satisfy /^(c) = | r'^{z)dz. The properties of In are well known. We only 
list some facts we are going to use. 

Lemma, (i) The following identity holds: 

{n + 6)/„+3 - 2uj{n + 3)/„+i - 2nc/„_i = 0, n G Z 

which implies 

T 2c ^ 80;^ ^ Suc^ {2c 32a;2x 
h = uh, I^ = -Io + —l2, ^6 = ^/0+ (y + ^)/2. (5.3) 

(ii) The integrals Jq and I2 satisfy the system 

AcI'q + 2w/^ = 3/0, 

4a;c/^ + (12c + 8a;2)/^ = 15/2. 

(iii) The ratio R{c) = /2(c) //q(c) satisfies the Riccati equation and related system 

c = 8c^ ~1~ AiUj'^c 

{8c^ + 4:UJ^c)R'{c) = -2ujc + AcR{c) + ujR^{c), . ' (5.4) 

R = -2ujc + AcR + uR^, 

which imply estimates 

2c , ^ c 3u , ^ 
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The equations in (i)-(iii) are derived in a standard way, see [9] for more details. 
The estimates flS.Sp follow from the fact that, in the (c, -R)-plane, the graph of R{c) 
coincides with the concave separatrix trajectory of the system (I5.4p contained in the 
triangle with vertices (0,0), {—^ijj^,u) and (0, |a;) and connecting the first two of 
them. Note that -R'(c) = — l/(2a;) at c = —\uj'^- We also use analyticity of the ratio 
R{c) at this point and properties of the phase portrait of (15.41) to verify the above 
statements. 

After this preparation, we turn to prove the estimates we used in the preceding 
sections. 

A2. The estimate for A2 = -Ijy^. By (g^]), (O and the first inequality in f l53|) . 
we have 

B2. The estimate for B2 = As /' = + (— + I' we have by 

llr'ip 6 15 \^ 5 15 y ^ ^ 



13. 4p as above 



_ 2 ucI'q + (3c + 3u;2) _ 2 wc + (3c + 2uj^)R 4w 
~ 5 2c/^ + Lol^ ~ 5 2c + ujR " T' 
To obtain the last inequality, we used both estimates in (15. 5p . 

C2. The estimate for C2 = X2 + Ao(l + ^) J'ffp . By ^12^ and (137[3|) we have 



Next, 



dj V {r^^o)"^ 

{r, V'o) = / [3r=^ - (3w + \o)r]dt = 6/^ - {6uj + 2\^)I[ = -2\oI[, 
Jo 

||r||2||^o||2 = J^r^dtJ^{3r^-3uj-Xofdt 

= 4/^ [9/^ - (18u; + 6Ao)/^ + {3u + Aq)'/^] 
= 4J^[(6c + (3w + Ao)2)/^ - (Gu + 6Ao)/^]. 

By (13. 5p . we have 

1 pT rK{k) 

I'^[c) = - r^dt = (rl/a) / dn'^{t)dt = V2roE{k). (5.6) 

2 Jo Jo 
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Making use of the identity E{k) = |7rF(|, — |, 1, k"^) where F is the Gauss hyperge- 
ometric function, we obtain an appropriate expansion to estimate E from above 



^ ' ~ 2 V 4 64 512 •••7 ' - - 4 V' 2 
with all removed terms negative. As I[ = -^^r, by (13. 5p this implies 

•2 lOwrg + 
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r/2 

-'2 — 



-^1 



/2 



Together with /^/^ > I'/, this yields 



kIPIIV'ol 



1 < 



6c + (Sou + Xof + —^{u + Ao)(w^ - lOurl + r, 



^0 



A2 f X2 + UJ 
^0 



4r2 



1 < 



A_2 / 73 

An 



where the equality is obtained by direct calculations. Therefore, 



2^2 



- 1 



with C2 > an absolute constant when d > 4 is fixed. 

As a by-product of our calculations, we easily obtain also the estimate 



An 



AqIIV^oII 



> 



X,(4 



1) + An 



D2. The estimates for D 



21 



and D22 = Mil • We proceed as in case Dl 



above. Making use of statements (i) and (ii) of the Lemma, we have 



^21 



_ _ 2^^4 ^ ^6)^^ ^ ^2 J/ _ 2ujr^ + 

-Acul'o + (18c + 7u^)l!^ _ -uIq + 



Similarly, 



^22 



5(2cJ^ + a;/^) 
2c/^ + w/2 _ (2c + u^)Io ^ 2c + 



< 5uj. 



I'2 



3/2 5/2 — a; Jo 3a; 

As before, we used that /2(c) //o(c) is a decreasing function and calculated the value 
/2(0)//o(0) = fa;. 
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